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Abstract. Bifurcation with spherical symmetry and mode interaction is known to 
produce robust heteroclinic cycles between axisymmetric steady-states. The case of 
mode interaction with representations of SO(3) of degrees £ = 1 and 2 has been 
studied in the context of the onset of convection in a spherical shell and central 
gravity force by [3] and [6] who explained certain numerical observations made in 
1986 by Friedrich and Haken [11]. This problem has recently regained interest thanks 
to the launch of an experiment of convection between two concentric spheres in the 
International Space Station, thus providing a system with (nearly) spherical symmetry 
(GeoFlow project, see [9]). In this case however, the onset of convection cannot be 
excited by modes of degree £ smaller than 3 due to technical constraints. Motivated 
initally by this experiment, we have analyzed the occurence of robust heteroclinic 
cycle in the case of £ = 3,4 mode interaction. This case is highly complex but, 
applying the methods of equivariant bifurcation theory, we have shown the existence 
of (generalized) robust heteroclinic cycles involving, not only axisymmetric states, but 
also and principally states with cubic symmetry. These objects are observable in the 
numerical simulations of the dynamics on the center manifold. We provide video clips 
showing the corresponding evolution of the pattern on the sphere. 



AMS classification scheme numbers: 34C37, 37C80, 37G40, 37L20 



1. Introduction 



Evolution equations which are invariant under a linear action of the group 0(3) are 
known to possess heteroclinic cycles, i.e invariant sets consisting of equilibria Pi, • • ■ , 
such that the unstable manifold of Pj is included in the stable manifold of Pj+i (mod 
k) in a structurally robust way (robust heteroclinic cycles are defined and discussed in 
[8]). The existence of these objects can be proved by bifurcation methods under certain 
conditions: 1) steady-state bifurcation occurs via mode interaction between irreducible 
representations of 0(3) of degrees i and i + 2) the coefficients of the quadratic terms 
of the equation on the center manifold satisfy certain relations which we shall make 
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precise later. It turns out that these two conditions are generically fulfilled when the 
initial problem is the (classical) model of Rayleigh-Bnard convection in a spherical shell 
(see [2]). The case when i = 1, which is the simplest situation, has been analyzed in 
details and has been fully solved in [6j. It was shown that heterochnic cycles in the above 
sense do indeed exist but are part of larger invariant sets which are called generalized 
heteroclinic cycles and which have been shown to be asymptotically stable. It was also 
shown that a robust heteroclinic cycle did still exist when the spherical shell was allowed 
to rotate (slowly) around an axis. 

Recently some interest has grown about the case of convection of a fiuid in a spherical 
shell with aspect ratio such that the onset of convection arises with mode interaction 
i = 3 and 4. This was in relation with the preparation of an experiment supported by the 
European Space Agency and named GeoFlow which consists of a spherical vessel filled 
with a fiuid subjected to an electrophoretic central force [9j. This device is placed in the 
International Space Station in order to simulate a self-gravitating fiuid. For technical 
reasons the aspect ratio of the spherical shell cannot be too small and the first unstable 
spherical modes are expected to be of degree £ = 3 or 4. Numerical simulations near 
onset, using a center manifold approach in the case of mode interaction, have shown 
fiow regimes which are quite suggestive of the presence of an attracting heteroclinic cycle 
involving steady states with cubic and/or tetrahedral isotropics [1]. Direct simulations 
show the same behaviour for some time but the fiow always stabilizes on a steady-state 
with tetrahedral symmetry after a few switches [13] . We interpret this contradiction by 
the fact that the equations for the modes with £ = 4 on the center manifold undergo 
a transcritical bifurcation of unstable steady-states (this is a well-known fact, see [7]), 
but have a quadratic term which is not very large and allows for a "bending back" of 
this branch with stable solutions above the turning point [12]. Then the equations for 
the modes with i = 3 generate the heteroclinic behaviour for the i = 4 steady-states. 
However it is likely that this "bending back" goes beyond the validity domain of the 
center manifold approximation, which would explain why the intermittent behaviour 
does not persist in the direct numerical model. In the Rayleigh-Benard convection 
model, this quadratic term has a coefficient usually very close to due to similar form 
of the gravity force and buoyancy force (see [2], [E]). This is a general fact valid for any 
even value of i and it has allowed for a description of a large family of heteroclinic cycles 
in the mode interaction for this problem [5]. However in the experimental device, 

the gravity field is replaced by an electrophoretic field which has a radial dependence 
in instead of for the buoyancy field, and this prevents the quadratic term in 
the i = 4 equations from being small (except possibly for special values of the Prandtl 
number). Moreover the heteroclinic cycles found in [5] are most probably always parts 
of more complex invariant sets and no stability assessment has been made except in the 
case i = 1 where it was precisely shown that the heteroclinic cycle is in fact embedded in 
a "generalized heteroclinic cycle" as was said above. To be complete it must also be said 
that in the classification of [5] (Theorem 5.5), the case i = 3 contains a mistake: indeed 
in this case precisely, a topological obstruction exists in the 2-dimensional space D4 © '^2^ 
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(set of points which are fixed by the subgroup 1)4® 1^2'^), preventing the existence of 
connections between axisymmetric states in this plane. This is one more reason to have 
a closer look at this case. 

Our aim in this paper is to investigate the occurrence of generalized heteroclinic 
cycles in the £ = 3,4 mode interaction under conditions which are compatible with 
onset of Rayleigh-Benard convection. When numerical simulations of the fiow on the 
center manifold are performed, a great variety of complex spatiotemporal regimes can 
be observed near bifurcation, depending on coefficient and parameter values. We have 
focused on the simplest possible situations. We show that robust generalized heteroclinic 
cycles connecting axisymmetric steady-states with ones having cubic symmetry can 
bifurcate. In [5], such objects involving steady-states with cubic symmetry where found 
to exist only for high values of t. 8, 9 and 13. Moreover these generalized heteroclinic 
cycles can be attracting. 

We begin by a description of the geometry of the 0(3) group action in the 3, 4 interaction 
and we describe the equivariant structure of the equations on the center manifold 
(coefficients listed in the appendix). Then we look at the primary and (some of the) 
secondary bifurcations of steady-states. In Section 3 we give a precise definition of the 
generalized heteroclinic cycles and we show the existence of such objects under certain 
conditions which are satisfied for an open set of parameter values. 

2. Bifurcation with £ = 3,4 mode interaction with 0(3) symmetry 
2.1. Isotropy types, equations on the center manifold 

A very basic step when studying a bifurcation problem with symmetry, is to determine 
the isotropy types of the group action on the center manifold and the corresponding 
fixed-point subspaces. Let us recall that given x G X vector space with a representation 
g — > T{g) of a group G {G = 0(3) in our case), we define the isotropy subgroup of x 
to be the subgroup {g E G / T{g)x = x}. The fixed-point subspace associated with an 
isotropy subgroup H is the set {?/ G X / T{H)y = y}. The fundamental property of G- 
equivariant evolution equations (or dynamical systems) is that given an initial condition 
Xq G Fix{H), the trajectory {x{t)} lies entirely in Fix{H). It is clear that if H G H' 
are two isotropy subgroups, then Fix{H) D Fix{H'). Also, two subspaces Fix{H) and 
Fix{H') are symmetric to each other, i.e. there exists an element g G G such that 
Fix{H') = T{g)Fix{H), if and only if H and H' are conjugated by an element in G 
(namely by g). The conjugacy class of an isotropy subgroup is called an isotropy type. 
From the above comments, it is clear that isotropy types are partially ordered. 

Natural irreducible representations of 0(3) of degree £ are spanned by the spherical 
harmonics ^£^"(^9, if) {—i<m< Therefore the dimension of the center manifold in 
the 3, 4 mode interaction is 16. Moreover, we know that the antipodal symmetry S acts 
trivially when i is even and as —Id when i is odd. We denote by the 2-element group 
of S. This is the notation which was introduced by Golubitsky, Stewart and Schaeffer 
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in [H] and which we will follow more generally to name the isotropy subgroups (see a 
correspondence table with Schoenflies notation in [8]). 

The diagram in Figure [T] below shows the tree of isotropy types for the £ = 3,4 mode 
interaction with 0(3) symmetry. This diagram has been obtained from the informations 
contained in [7j where notations of the groups are also explained. The numbers in 
parenthesis on the left side indicate the dimensions of the fixed-point subspaces on each 
corresponding line. More on the structure of these subspaces will be given in next 
sections. 

The isotropy types with one dimensional fixed point subspaces are the "pure" i = 4 
mode isotropics 0(2) © (axisymmetric modes) and 0©Z2'^ (cubic symmetry modes). 
We know from the Equivariant Branching Lemma [H] that corresponding branches of 
equilibria bifurcate from the trivial state. In the "generic" case, it is known that the 
branches are transcritical (presence of quadratic terms in the bifurcation equations) and 
unstable [7]. We shall be interested in the case when the pure quadratic terms have a 
vanishing or nearly vanishing coefficient, a situation which occurs in classical Rayleigh- 
Benard convection. In this case, not only these bifurcated branches can bend back and 
gain stability, but also other types of solutions with lower isotropy can bifurcate through 
secondary and tertiary bifurcations. The unfolding of this quadratic degeneracy has 
been studied in details in [12]. The complete picture looks already quite complicated 
in the pure i = 4 case. The coupling with i = 3 modes makes the problem even 
more complicated, as can be already seen from the diagram of isotropy types. We 
shall therefore focus on some specific, relatively simple situations where the existence 
of heteroclinic orbits can be proved and where the number of steady-states involved in 
the construction of heteroclinic cycles is limited. 

Let X = (a;_3, ■ ■ ■ , X3) denote the component along the i = 3 spherical harmonics 
, and Y = ■ ■ ■ ^y^) denote the component along the £ = 4 spherical harmonics 
Y^ . Note that, Y^"^ = (—1)™!^™' and therefore the same holds for the coordinates (real 
space). We list in the table below the components of the fixed-point subspaces with 
dimension lower than 6. These isotropy types are the important ones for the subsequent 
study. The reader can easily complete the list by looking at the tables in |7J. Our 
convention is that coordinates which are not indicated should be set equal to 0. Also, 
we write (general notation): z = Zr + izi E C. 

Let us now write the equations on the center manifold. Only the leading terms 
need be considered here. The context of this paper is the onset of convection in a self- 
gravitating spherical shell with uniformly heated inner boundary. This assumption, as 
shown in [5] (in [2] in a more general context), has strong consequences on the signs 
of coefficients of the leading part of the equations on the center manifold. Because S 
acts non trivially when i = 3 but trivially when i = 4, the "pure" mode, 9-dimensional 
subspace associated with £ = 4 is fiow- invariant (this is Fix{'L2^))- We perform a generic 
change of variables in the parameters if necessary, so that there are two free parameters 
Hi and /i2 such that the C = 3 modes become unstable when /ii > and the £ = 4 modes 
become unstable when /i2 > 0. 
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Figure 1. The isotropy types for the £ = 3,4 mode interaction. Numbers on the left 
indicate dimension of corresponding fixed point subspaces 

Then the leading part of the equations in these coordinates read 

Xj = niXj + l3Qf\x-3, ■■■,X3, y-4, • • • , ^4) + Till-'S^ll^ajj + 72Cj^^(a;_3, • • • , ^3) (1) 
ilk = IJ'2yk + bQf^{x-3, • • • , ars) + cQ^^\y-A, • • • , ^4) + di WW^Vk 

+ d2Cf (y_4,---,2/4) (2) 

where the Q's denote quadratic terms and the C's denote cubic terms which are given 
exphcitly in the annex. In the context of Rayleigh-Benard convection in a spherical 



Robust heteroclinic cycles in spherically invariant systems 



6 





Fir(H) 


l^o\J Vi L KJ VJ y ±± 


J. t/Jb I ^ 




yo 


U2 to <2-'2 


\ 1 lr\ 1 In lit 1 

I,y0) y2r; y4r j 




7//1 = ~l~71?/n 
y4r -'-'lyi} 




V-^iJ; yO) y4ry 




i ^ r-. 1 fr\\ 


1^3 


1,J'3«5 yo; y3r J 


D4 © 


(2/0,l/4r) 


D2 


(a;2i, yo, y2r, y^r) 




(a;3r,2/o) 


Di 


{Xo,X3r,yo,y3r) 


030^2" 


(2/0,l/3r) 




{Xo,X2r,yo,y2r,yAr) 


0- 


(a;2j, l/4r = Wo) 


Z2 © Z2" 


(1/0, Z/2, 2/2, 1/4, 2/4) 


D| 


{xo,yo,y4r) 


Z4" 


(x2,X2,2/0,2/4,2/4) 



Table 1. Fixed-point subspaces for representatives of the isotropy types such that 
dim Fix{H) < 5. We have noted 77 = Another useful representation of 

FixiO © ^2") is ysr = i^yo with v = J^. 



shell, the coefficients satisfy certain relations which, using the normalization applied to 
compute the equivariant polynomials listed in the Appendix, are as follows 

Hypotheses (HI): 

(i) h ~ —2/3. We do not loose generality by assuming that /3 = — 1 (after a suitable 
change of time scale). 

(ii) 71 + 2O72 < 0, 

(iii) + ^2 < 0, j = 1,- ■ - ,4, 

(iv) c ~ will be considered as an additional parameter. In all the following we restrict 
to the case c > 0. 

2.2. Primary branches of steady- states 

There are two primary bifurcated branches of steady-states which correspond to the 
maximal isotropy types (those with dimFix(if) = 1). We detail these solutions and 
compute the eigenvalues of the Jacobian of the vector field at these solutions, in order 
to study later their stable and unstable manifolds. 

2.2.1. The axisymmetric equilibria On the axis Fza;(0(2) © ^2^) the (scalar) 
bifurcation equation reads 

= ^^2yo + 9c?/o + {di + ^2)2/0- 

It follows that the bifurcated equilibria satisfy the relation 

/i2 = -9c2/o - {di + d2)yl (3) 

A turning point exists at 7/0 = ~ 2{(h+d2) ^^'^ Il-(iii) the parabola is oriented towards 
2/0 > in the (/i2, 2/o) plane. Even in the hmit c = 0, the two branches of the parabola 
correspond to symmetrically distinct states. 
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We denote by a± these branches of axisymmetric equihbria along the 
component. They generate two-dimensional 0(3) orbits of equilibria. In order to 
determine their stable and unstable manifolds we need to compute the eigenvalues of 
the Jacobian matrix L^^ of the vector field linearized at a±, hence to determine first 
the isotypic decomposition of the action of 0(2) © ^2^^ in the 3, 4 representation of 0(3). 
A straightforward (and classical) analysis shows that L^^ decomposes into two 1x1 
and seven 2x2 diagonal blocks along the coordinates Xq, yo, {xj,Xj) (j = 1,2,3) and 
{Vm, Vm) ijn = 1, 2, 3, 4) respectively. The block in the {yi, yi) subspace is the matrix, 
corresponding to the fact that this plane is tangent to the 0(3) orbit of a± [7J. The 
other eigenvalues are listed in the following table [2] where yo is the coordinate of a±. 
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Table 2. Eigenvalues of the linearized vector field at q = a± 



2.2.2. The cubic equilibria Let us consider the axis Fix{0®1>2^) defined by the relation 
?/4j, = A/5/14t/o- The bifurcation equation along this axis parametrized by yo reads 

2 /12 , 16 , , o 
= /i2l/o + 14c?/o + (ycti - ^"2)2/0 

and therefore that bifurcated branches are given by 

/i2 = -14cyo - (yc^i - ■^d2)yl (4) 
We shall assume: 

Hypothesis (H2): 3di - 1^2 < 0. 

(H2) combined with (Hl-iii) is equivalent to rfi < and 21/4(ii < ^2 < — rfi. Then the 
parabola is oriented towards > in the (/i2;2/o) plane. Even in the limit c = the 
two branches of the parabola correspond to symmetrically distinct states. 

We note these solutions /3 = /3+ or They generate three-dimensional 0(3) orbits 
of equilibria. The eigenvalues of the Jacobian matrix Lp^ of the vector field linearized 
at (3± are listed in the following table [3l together with the corresponding eigenspaces. 
Here again the multiplicity and eigenspaces follow from the isotypic decomposition of 
the action of O © 1^2^ in the £ = 3,4 representation of 0(3). The convention is that 
coordinates which do not appear in the definition of eigenspaces must be set equal to 0. 
The j3± solutions are parametrized by their coordinate yo- 
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Eigenvalues 


Multiplicity 


Eigenspaces 


(t'o = Hi + 12yo 
o-f = A^i - 6yo 


1 
6 


r ^ vi5 ™ T 
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3 
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A^ = -40cyo + l0/7d2yl 


3 


{y2i,y3 = \/72/i} 



Table 3. Eigenvalues of the linearized vector field at P = P± 



2.3. Suhmaximal bifurcations 

There are two kinds of branches of steady-states with non maximal isotropy in this 
problem: 1) secondary branching in the pure £ = 4 mode subspace, due to the quadractic 
degeneracy c ~ 0; 2) branching occurring in mixed mode subspaces (with /ii as a natural 
bifurcation parameter). The first case, which is independent of the mode interaction, 
has been studied in detail by [12]. Equilibria with Dfc © ^2'^, A; = 2, 3, 4, have been shown 
to bifurcate in certain regions of parameter space. This combined with the bifurcations 
in the £ = 3 directions makes a very complicated bifurcation diagram which we do not 
intend to describe here. Our aim here is rather to focus on the simple case of suhmaximal 
bifurcation, i.e. branching of solutions with suhmaximal isotropy [8], [7j. In our case, 
this means branching in one of the invariant planes listed in Table [H 



2.3.1. Bifurcations in the pure i = 4 invariant planes According to Fig. [T] there are 
two kinds of such planes. Let us choose the representatives given in Table [H i.e. with 
coordinates resp. Fix(D4©Z2'^) = {yo,l/4r} and Fix(D3©Z2'^) = {yo,y3r}- Each of 
them contains the axis Fix{0{2) 0)1^2^) = {Vo} and two copies of the axes of cubic 
isotropy O © {?/4r = ±mo} in Fia;(D4 © ^2") and {ysr = ±Jyyo} in FixiBs © Z2') 
(where rj = \J^, v = \J^)- 

It is a known fact that, when the quadratic coefficient c is fixed non zero, the 
bifurcated branches with maximal isotropy (types a and (3) are transcritical and no 
other solutions occur [7]. However we assume the degeneracy hypothesis (Hl-iv). Let 
us focus first on the plane P = {yo,y4r}- Writing u = yo and v = y^r to simplify 
notations, the equations read 

26 

u = ^^u + c(9m^ + lAv'^) + diu{u^ + 2v'^) + d2u{u^ - —v"^) (5) 

13 30 

V = ^i2V + lAcuv + div{u^ + 2v'^) + d2v{ — + —v"^) (6) 

We note (3± (resp. (3±) the equilibria on the invariant axis {y^r = VVo} (resp. 
{Uir = ~Vyo})- Equilibria off the invariant axes (hence with D4©Z2'^ isotropy) must 
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satisfy the conditions 
7c_ 
4do 



u 



1^2= - 



and 

(7^1 + 43^2) - (2c/i + —d2)v^ 

2 ' 



(7) 
(8) 



We see that the bifurcation from a (axisymmetric) equihbria in the plane P is a pitchfork 
and occurs at the /i2 value 



7c^ 
16^2' 



When c > and (i2 < the branching occurs off the q;_ equilibrium. Let us denote by 
7, 7 these equilibria. As /i2 is increased, they move closer to the invariant axes with 
cubic isotropy until they cross them at the value 



//f = -(21rfi + 94^2) 



This corresponds to a transcritical bifurcation from j3 and j3 solutions. When c > 
and c?2 < 0, the secondary bifurcations occur from the «_ branch (at /X2) and from the 
branches /3_ and j3- (at /if). The latter correspond to the vanishing of the eigenvalue 
A2" (see Table [3]). 

The phase diagram in the plane P when < ^2 < fJ'2 looks like Fig. [21 





/ ^ 






+ 


p-l // 





n 



Figure 2. Phase portrait in P ~ Fia;(D4 ® 22^^) when < ^2 < ^2- Equilibria /3± 
and /3± are exchanged by rotation (/) ~ i: /A. 



An examination of the eigenvalues in Tables [2] and [3] show that as long as /i < /Xg, 
there are no other secondary bifurcations in the pure £ = 4 modes. It can be shown that 
indeed no other branch of equilibria exist in this case [I2]. This means in particular 
that the secondary bifurcations with D3 © ^2^ isotropy, which can be described in a very 
similar manner in the plane (|/o,l/3r), occur at values of /i2 larger than /if- 
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More precisely, let us now consider the dynamics in Fix(j)3 © '^2'^)- We use the notation 
yo = u,y:ir = w. Then, 

17 

ii = H2U + c{9u^ - 21w'^) + diu{u^ + 2w^) + d2u{u'^ - j^^^) (9) 

17 5 

w = iJ,2W — 21cuw + diw{u^ + 2w^) + d2w{ w"^) (10) 

28 56 

The invariant axes with cubic isotropy are defined hj w = ±^/^m and the branches of 
equilibria on these axes is given by 

/27 , 36 , , 2 9,12 , 16 , , 2 

^i2 = 2lcu - [—di - —d2)u = 2lcu - -{—di - —d2)u (11) 

The bifurcation scenario is similar to the one in the plane P, except for the following. 
Comparing the above expression with (jl]), we deduce that the solutions with m > (resp. 
■u < 0) belong to the orbit of /?_ (resp. /3+). Moreover, the secondary bifurcations occur 
when the eigenvalues A3 or A3 vanish, and this happens at values of /i2 which are larger 
than /i2- In particular, when is smaller than 

/i^ = -8(168cii + 17ci2);T2 
there are no secondary equilibria and the phase diagram looks like in Figure [3l 



w 




Figure 3. Phase portrait in Fia;(D3®Z2°) when ^2 < A"- Notations /3± mean that 
the equilibria belong to the group orbits of j3±. 



2.3.2. Invariant "mixed-mode" planes According to the isotropy tree in Fig. [T] there 
are three such planes: Fix{0{2)~), Fzx(Dg) and Fix{Q~). The first two planes contain 
the invariant axis Fix{0{2) © '^2^) while the third plane contains the axis Fix(0©Z2'^). 
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Using the expressions given in the Appendix for the quadratic equivariants, it is easy 
to check that the conditions of apphcation of Proposition 1 of [2] are met when (HI) is 
satisfied. Therefore heterochnic orbits connecting the two bifurcated equihbria lying on 
the invariant axis exist for an open set of parameter values. 

Let us be more specific about the plane 11 = Fix{0~), which can be defined by the 
coordinates X2i and y^r = VUo iv = ■\/5/14). The axis of symmetry {y^r = VUo} contains 
the equilibria Using j/g as a parameter for this axis, the equations for the flow on 11 
read 

X2i = f^iX2i + l2x2iyo + (27i + 2072)a;^i (12) 

m = /i2l/o + Ucyl - + (12/7rfi - 16/49rf2)z/o (13) 

The bifurcation scenario for this kind of planar system is classical and we don't go into 
details of the calculations (see [3J). Let us assume first that c = and fix /i2 > in 
order to insure existence of branches with maximal isotropy (thanks to HI and H2). We 
set 12/7cii - 16/49(i2 = D and 

[D < hj hypothesis H2). Then as fii is varied the following events occur: 

- For < /ii, there exist no equilibria off the invariant axis and (3± are both sinks in 
H. 

- At /ii = /ii a pitchfork bifurcation occurs from /?+ off the axis. We note S and S the 
bifurcated equilibria with isotropy 0~ (which are exchanged by taking the coordinate 
X2i to its opposite). 

- As /ii is increased further, a Hopf bifurcation occurs from 6 (and 6). 

- The two limit cycles grow until they collide with the stable manifold of the origin. 
Then they disappear and two saddle-sink connections from /5+ to /3_ are established at 
a value /ii = fii which satisfies /ii < /ti < 0. 

- The equilibria 6, 6 die off at the origin in a "reverse" bifurcation when fii = 0. 

- The heteroclinic orbit from /5+ to /3_ persists until a new pitchfork bifurcation occurs 
from P_ at a positive value of fii. 

These properties persist when c is taken non zero but small. 

Similar scenarios occur in Fix{0{2)~) and in Fixl'Of), with secondary bifurcation 
from the axisymmetric states a±. It can be easily checked that heteroclinic orbits 
between these states can coexist in both planes, however the eigenvalues at a_, resp. 
a+ along these connections, which are given by Aq and A3, have the same sign, and 
therefore these orbits cannot form heteroclinic cycles. In the next section we shall 
consider situations where such connections do not exist. 

3. Robust generalized heteroclinic cycles 

As it was already shown in [6j, the usual definition of a heteroclinic cycle (as given 
in [8]) is not sufficient to describe recurrent heteroclinic behaviour in systems with 
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0(3) symmetry. Indeed it may happen that instead of a cycle, the set of heterochnic 
orbits makes a non trivial (closed) network with multidimensional connections between 
equilibria. Hereafter we give a suitable definition of a robust generalized heteroclinic 
cycle in order to take account of this generality. Given an equilibrium ^ of ([I])-®, we 
denote by its 0(3) orbit and more generally, given a submanifold M in phase space, 
we set OM = {g ■ M \ g G 0(3)} where g ■ denotes the action oi g G 0(3). We 
also denote by W^^ (resp. W^) its stable (resp. unstable) manifold (see [8]). Note that, 
OVF|'" = Wq^. The following definition is stated for 0(3) equivariant systems but is 
obviously generalized to any other group of symmetry. 

Definition 1 A set of 0(3) orbits of saddle equilibria 0^i,...,0^„, defines a 
generalized, 0(3) invariant heteroclinic cycle (a GHCj if for any j G {1, . . . , n}, 

The GHC itself is the union IJ^- Wq^. ■ This definition is stronger than others because we 
request that all unstable trajectories be heteroclinic orbits in the GHC. On the other 
hand the manifolds of connections can be multidimensional and can follow multiple 
circuits. 

We look for heteroclinic cycles which are robust under small 0(3)-equivariant 
perturbations (equivariant structural stability). This can indeed happen if the following 
is true. 

Condition for structural stability: whenever two equilibria ^ and ^' in the GHC 

are connected, i.e. satisfy fl W^, ^ 0, then there exists an isotropy subgroup H such 
that the following properties hold (transversality conditions): 

(i) i,i'eFix{H) 

(ii) either ^ is a source or ^' is a sink in Fix{H), 

(iii) n W^, n Fix{H) ^ 

Our aim in this section is at finding such objects. We begin with the case when only 
a and f3 equilibria are involved. We look first at heteroclinic connections in Fix{'L2^) 
(pure £ = 4 modes), then at connections involving both £ = 3 and 4 modes (so- 
called "mixed modes"). Remark that our study does not exhaust all possible types 
of generalized heteroclinic cycles. In fact it is easy to find other types which can be 
deduced from those analyzed here and which however do not add better understanding 
of the problem. 

3.1. Generalized heteroclinic cycles with a and P equilibria 
3.1.1. Connections in Fix{Z2^) 

Proposition 1 Suppose hypotheses (HI) and (H2) are satisfied and in addition < 
ji2 < IJ'2- Then 0^_^ (the 0(3) orbit of equilibria P+) is attracting in Fix{Ij2^) and 
moreover the equilibria a± and (3± satisfy the properties (i)-(iii) above. 
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Proof. We rely on [12] to assert that under the conditions of the proposition, there are 
no secondary branches of equihbria. Moreover, the system is determined at cubic order, 
and this truncated vector field is a gradient (see [7j). Therefore there are no other non 
trivial bounded sets than a and (3 equilibria and heteroclinic orbits connecting them 
in Fix{Z2^) (note however that this gradient property prevents heteroclinic orbits to 
form a cycle.) Now, examination of the eigenvalues A" and (tables [2] and [3]) show 
that under the conditions (HI) and (H2), /5+ is stable as a relative equilibrium (i.e. 
its 0(3) orbit is stable) and there are no other stable equilibria. The second part of 
the proposition follows from the fact that the dynamics is bounded in Fix{Z2^) when 
conditions (H1)-(H2) are satisfied. This in turn follows from the fact that the quadratic 
form + ■ C^'^\Y) is negative, where C^'^\Y) is the cubic map defined in 

Table lA5l The proof is straightforward if one assumes \d2\ small. In fact, thanks to the 
assumption that c is close enough to 0, we can apply a theorem by M. Field which 
states the existence of a flow-invariant and attracting topological sphere in Fix{'L2'^)- I 
Let us now be more precise on the phase portrait. Restricting to the planes 
-Fza;(D4 © ^2^^) and Fix(D3 © ^2^^), the pictures are like in Figures [2] and [31 It is 
more delicate to describe the phase portrait in higher dimensional fixed-point subspaces 
but it is instructive to consider the 3 dimensional Fix(D2 © ^2^^)- This space contains 
three copies P, P' and P" of Pix(D4 © ^2*^) which intersect all three on the same axis 
L = Fix{0 Q) 1^2^)- Using coordinates {yo,y2r,y4r) (see Table [I]), these planes have 
equations y2r = and y2r = M^Vo - ^Vir) and L = {?/2r = 0,y4r = V'5/l%o}. 
The three planes are exchanged by the action of A^(D2)/ D2 — -D3 (we note iV(D2) the 
normalizer of D2 in SO (3), which is known to be O, see [H]). This structure is sketched 
in Fig. m below. Assuming the existence of the fiow-invariant sphere, the bounded 
dynamics restricts to that sphere and can therefore be easily seen on a 2D picture (see 
Fig. [5]). We see that /3+ and /3+ are sinks, /3_ and have robust connections to 
a_ (in Pix(D4 © ^12'^)) which itself has a robust connection to This combined with 
the phase portrait in FixiD^ © ^2^^) (Fig. [3]) gives a good description of the various 
connections in Fix{'L2'^)- 




Figure 4. The subspace Fix{'D2®'^2')- 
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3.1.2. Connections involving £ = 3 modes Since /3+ is a sink in Fix{'L2^), we look 
for its unstable manifold involving necessarily £ = 3 perturbations and such that 
the conditions (i)-(iii) are met. We already know from Section 12.3.21 that saddle-sink 
connections between j3+ and (3- exist in the plane 11 = Fix{0~) when /ii > fii where 
fii < 0. However there are constraints in order to insure that conditions (i)-(iii) above 
are satisfied. If yUi > 0, the eigenvalues a"" with j = 1,3 are positive as well as af", 
and then the problem becomes too complicated (although simple attractors like relative 
equilibria or relative periodic orbits can be numerically observed at certain parameter 
values). We therefore suppose /ii < 0. Then clearly, from the expressions in Tables [2] 
and [3l the following is true: aj"*" < ( j = 0,1,3), < 0, a^' < 0, cxf''' < and 
a?- < 0. 

The following lemma shows that an open interval exists in which all the other eigenvalues 
are negative except cTq^ > which is a necessary condition for the saddle-sink connection 
in n. Let and y^ denote the values of yo at the equilibria a±, resp. P±, in Tables [2] 
and [31 

Lemma 1 Assume (H1)-(H2) and /ii < 0. Define fi[ = —7y°, fJ^" = —'2y^ (negative 
values). The following estimates hold: 

(i) if < min(/i']^, /i'/), then aj~ < ( j = 0,1,3), < 0, erf" < and a2~^ < 0; 

(ii) if Hi > 6/i'/, then ctq^ > 0. 

(Hi) There exists a value ^2 < such that the following is true: if21/Adi < d2 < d2, 
then /i" < fi[ while if > d2 > d2, then Gfi" < fi[ < fi". 

Proof. By Tables |2] and [3], the signs of eigenvalues in point (i) are satisfied if 
/ii < Qy- ,fii < -7y+, /ii < Qy-, /^i < -2?/J 



and point (ii) is satisfied if 

/ii > -12y: 
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The point (ii) is therefore immediate. Solving (jl]) and for y^, we find that 



(14) 



y± = ^ (-7c T V49c2 - Dfx^) 



(15) 



where D = 12/7(ii — 16/49^2- Recall that by our hypotheses, D < and di + d2 < 0. 
Replacing these formulas in the above inequalities and using our hypotheses, point (i) 
follows without difficulty. To prove point (iii) one can also easily check the following: if 
d2 = then G/i" < fi[ < while when d2 21/Adi (from above), then //'^ is bounded 



The value corresponds to fii in Section [2.3.21 The determination of the lowest value 
fii at which the saddle-sink connection (3^ — > /?_ is established, is very cumbersome. 
However in all the tests which have been done, it is observed numerically that fii < /i'/. 
This leads us to making the following ansatz: 

Ansatz 1 In the hypotheses of LemmaUi Ai < f^'i- 

It follows from lemma 1 that when ^2 is close enough to 0, fi[ < fi'l, therefore when 
fii lies between these two values the eigenvalue cxg^ is positive while the saddle-sink 
connections /5+ — /?_ exist in 11 and all the other eigenvalues in the i = 3 directions 
remain negative. It turns out that cxg''' corresponds to the eigendirection transverse to 
the plane P = Fix{D4 © '^2^) in the subspace Fix{Df) = P©n (see Tabled]). Therefore 
the unstable manifold of a+ in FixlDf) is two-dimensional, with eigendirections y^r 
(eigenvalue A4''') in P and X2i (eigenvalue ct2^) transversally to P. This situation is very 
similar to the one observed in the i = 1,2 mode interaction [6j. It can be shown with 
the same proof that the flow in the space FixlDf), for d2 and c in a neighborhood of 
0, looks like in Fig. [6l In fact, as in the Proposition 4.6 of [6], one can even show the 
following 

Proposition 2 In addition to hypotheses (HI) and (H2) suppose that d2 is close enough 
to and the Ansatz 1 is true. Then there exists an open interval of (negative) values of 
fii such that the full (4- dimensional) unstable manifold of a+ is included in the stable 
manifold of the group orbit of f3-^ in the 5-dimensional space FixCL^'^). 

Proof. We exploit the fact that if ^2 = 0, then the problem in Fix{Z4~) is identical to 
the problem restricted to Fix{'L2~) in the £ = 1, 2 interaction case which was studied in 
[6]. Indeed the geometry is similar, in particular we note that A^(Z4~)/ Z4" ~ 0(2) acts 
in the same manner in FixCL^") than N{'L2~)/ does in the £ = 1, 2 case. Moreover 
there is only one pure ^ = 2 cubic equivariant and if ^2 = the same is true in 

the £ = 4 case. The equations in Fix{'Li~) then read (at leading order): 



and fii 



—00. 



I 



X2 = fiiX2 + 7x2yo - V70x2yA + (271 + 2072)^2X2 

yo = /i22/o + c{9yo + Uyj) - 14x2X2 + diyo{yl + 2^42/4) 

2/4 = /i22/4 + 14cyo2/4 + -^xl + diyi{yl + 22/42/4) 



(16) 
(17) 



(18) 
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The proofs of Propositions 4.6 and 4.7 in [6j apply to this system, and taking ^2 close 
enough to does not modify the result. Numerical computations show evidence that 
the result persists for a larger range of values of d2 as well as of c. I 




Figure 6. The phase portrait in Fix{'Df) (in grey the heteroclinic connections in P). 



3.1.3. The generalized heteroclinic cycles We can summarize the previous results in 
the following theorem: 

Theorem 1 Assume hypotheses HI and H2 and the ansatzUl are true. Then, for open 
intervals of parameter values /ii < and fi2 > 0, a generalized heteroclinic cycle exists, 
which connects the 0{3)-orbits of equilibria of a and j3 types and whose structure de- 
pends on the values of the coefficient 

(i) There exists a neighborhood I of in M. such that when d2 & I , all four types a± and 
I3± are involved in a structure which is described in Propositions [I] and [H and sketched 
in Figure^ 

(ii) There exists a value d2, 21/Adi < d2 < 0, such that when 21/4(ii < ^2 < d2, the 
generalized heteroclinic cycle has a simpler structure than in (i), involving a+ and j3± 
only, according to the diagram of Figure 0. 

Proof of the Theorem. Point (i) is a direct consequence of Propositions [1] and 
El To prove point (ii) note first that equilibria 7 and 7 with isotropy D4 © ^2^^ bifurcate 
transcritically (and unstably) from /5_ when ^2 > ^2^ as shown in Section 12.3. 1[ Then 
the phase portrait Figure [5] (in Fia;(D2 © ^2^)) is replaced by Figure [9] below. The 
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Figure 7. Scheme of the heterochnic cycle when c?2 close to 0, showing the eigenvalues 
corresponding to the (manifolds of) connecting orbits 




Figure 8. Scheme of the heteroclinic cycle when d2 < d2 and /i2 > ^2 ' showing the 
eigenvalues corresponding to the (manifolds of) connecting orbits. 

heteroclinic cycle does not involve q;_ anymore. Moreover if d2 is large enough (in 
absolute value), an interval of values of /ii exists such that all eigenvalues involving 
i = 3 eigendirections are negative but ctq^ which is positive due to the existence of the 
heteroclinic orbits in 11. In this case the heteroclinic cycle has a much simpler structure. 

I 




Figure 9. Fia;(D2 © ^2*^): the phase portrait on the portion of invariant sphere 
bounded by P and P" (see Fig. [SJ in the case /i2 > ■ 
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3.2. A simple heteroclinic cycle involving mixed-mode equilibria 

We have seen in Section 12.3.21 that when /ii < /ti < 0, the presence of a stable fixed 
point 6 (and its symmetric 6) in the plane 11 prevented the existence of the saddle- 
sink connections in that plane. However if we now look at the 3-dimensional subspace 
Fixijyf) = P©n, a saddle-sink connection can exist in that subspace between 6 (in 11) 
and /3+ (in P). The process by which such a connection is created is the following: first, 
a steady-state bifurcation from 6 occurs off the plane 11 at a value /i'^ > fii. This new 
steady-state with D4 symmetry is a sink. As /ii increases, it travels toward the plane P 
until it disappears in Then the saddle-sink connection is established. A completely 
similar situation holds in the i = 1,2 interaction and was analyzed in more details in 
[3]. We refer to this paper for the formal proof of existence of this connection in the 
limits a = 0, d2 = 0. 

Once the saddle-sink heteroclinic connection exists, a heteroclinic cycle is set since 
/3+ belongs to the group orbit (9/5+ (see the diagram [19] below). Moreover this cycle is 
simple: it involves unstable manifold of dimension 1 only (in Fix{0~)). This was called 
a "type II" heteroclinic cycle by |3j. 

Of3+ 05 ► 0(5+ (19) 

Variants of this cycle involving periodic orbits in II do also exist for suitable 
parameter values. 

3.3. Stability and numerical simulation of the heteroclinic cycles 

The asymptotic stability of a " simple" heteroclinic cycle depends on the relative strength 
of the stable and unstable eigenvalues of the Jacobian matrix calculated at the equilibria 
which compose the cycle. Of particular relevance are the eigenvalues with eigendirections 
along the heteroclinic connections. General stability conditions for such cycles were 
derived by p^. These conditions apply to the simple cycle (fT9|) . but not to the 
generalized heteroclinic cycles of Section 13.1.31 Indeed one of the conditions is that 
for each equilibrium in the cycle, its unstable manifold consist of points with the same 
isotropy type. This is clearly not the case for the cycles in Figures [7] and [8j for example 
the unstable manifold of /3_ in the case of figure [H] intersects the strata of isotropy 
type D3 © '^2' (connections to equilibria of type Q!_|_), D2 © '^2^ and (connections to 
equilibria of type /5+). 

A similar situation holds in the £=1,2 mode interaction CclSG, clS shown in [6], however 
there the geometry was simpler and stability conditions could be derived by construction 
of a first return map using local orbit space reduction. 

Here we are facing two kinds of difficulties. First, the generalized cycles are in fact 
networks of connections between equilibria, in which different itineraries are possible due 
to the multidimensionality of the connecting sets. Hence the question of which itinerary 
is indeed followed by the asymptotic dynamics arises. This problem was investigated 
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by [T] in a simpler context and no general answer exists so far. Secondly, the unstable 
manifold associated with the same eigenvalue A3" (of multiplicity 3) is foliated with 
heteroclinic connections to a+ equilibria as well as to /5+ equilibria, as indicated on the 
diagrams of Figures [7] and [81 Therefore it is likely that nonlinearities play a role in the 
asymptotic dynamics near these heteroclinic cycles. 

The issue of finding the asymptotic stability conditions and minimal attracting sets 
for generalized heteroclinic cycles is therefore a difficult one and it goes beyond the aim of 
this paper, which is to present a new (and relevant) type of heteroclinic cycles occurring 
in systems with spherical symmetry. This issue will be the topic for a future work. 
Below we simply show on some examples the numerical evidence for asymptotic stability 
behaviour of the three types of generalized heteroclinic cycles which we have described 
in the previous section. These simulations have been realized on a Matlab code for the 
integration of the cubic order system. Videos showing the time evolution of the pattern 
associated with these time series can be seen on http: //www. math. unice . f r/~chossat[ 

3.3.1. Asymptotic dynamics near the heteroclinic cycle of Figure This cycle is 
numerically observed for a non empty set of values of the coefficients and parameters in 
the equations. Recall that it is assumed that ^2 (as well as c) is small in absolute value. 
In the following calculations we have taken: 

/5 = 0.08, 6 =-0.09, c = 0.01, 71 = -0.5, 72 = -0.5/36, rfi = -1, rfg = -0.008. 

The bifurcation parameter values are /xi = —0.09 and 112 = 0.03. Then a typical time 
series is shown in Fig. [101 This figure shows the time series (after transients have been 
eliminated) for the energy of the modes with £ = 3 (in blue) and 4 (in red). It can be 




observed on the video that the dynamics passes near the axisymmetric states of type 
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a_, but never close to those of type a+. When \d2\ is larger the flow between /3_ and 
/3+ follows the connection in D3 © '^2'^ and does not come close to axisymmetric states 
anymore. 

3.3.2. Asymptotic dynamics near the heteroclinic cycle of FigurelE We take the same 
coefficient values as in the previous subsection, except that now ^2 = —1-6. In this case 
the secondary steady-states 7 exist and the cycle of Figure [H] is realized. The bifurcation 
parameter values are fii = —0.1 and fi2 = 0.03. Then a typical time series is shown 
in Fig. [m On the corresponding video the axisymmetric steady-states are clearly not 

0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 


5000 6000 7000 8000 9000 10000 

Figure 11. Time series of energies associated with the modes £ = 3 (in blue) and 4 
(in red) in the case of Figure [51 

present. 

3.3.3. Asymptotic dynamics near the simple heteroclinic cycle of diagram / f75|) We take 
again d2 = 0.08, all other coefficients identical to the previous ones, but now fii = —0.13. 
In this case the /5+ —>■ /3_ connections in Fix{0~) are broken by mixed-mode steady- 
states which undergo saddle-sink connections in Fix{Df) to the /3+ equilibria. A typical 
time series is shown in Fig. [T2j 
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Figure 12. Time series of energies associated with the modes £ = 3 (in blue) and 4 
(in red) in the case of diagram (fTH)) . 
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Appendix A. The quadratic and cubic equivariant maps 

In the tables below we list the coefficients of quadratic and cubic maps in equations ([H) 
and ([2]). Conventions are as follows: Number m in the "Equation" column indicates the 
index of the component along the spherical harmonic Y^, with £ = 3 or 4 depending 
on the map under consideration. We list only the coefficients for the components with 
m = 0, ■ • ■ , +i, since the components for negative m are obtained from the ones with 
positive m according to the rule Yf'^ = {—1)"^Y^. The numbers in parenthesis in 
the column "terms" indicate the indices of the corresponding terms with the following 
convention: for the quadratic terms Xjj/j (—3 < z < 3, -4 < j < 4), [ij) for the 

terms XiXj (—3 < « < j < 3) or y^yj (—4 < « < j < 4), and (ijk) for the cubic terms 
ytyjyk (-4 < i < j < /c < 4). 



Equation 


Term 


Coefficient 


Equation 


Term 


Coefficient 


3 


(0,3) 


-3^/7 


2 


(-2,4) 


V70 




(-1,4) 


V42 




(-1,3) 


-VTi 




(1,2) 


3V6 




(0,2) 


-V3 




(2,1) 


-Vso 




(1,1) 


4V2 




(3,0) 


3 




(2,0) 


-7 










(3,-1) 


Vso 


1 


(-3,4) 


V42 





(-3,3) 






(-2,3) 






(-2,2) 


-V3 




(-1,2) 






(-1,1) 


-VTE 




(0,1) 


VT5 




(0,0) 


6 




(1,0) 


1 




(1,-1) 


-VT5 




(2,-1) 


-4%/2 




(2,-2) 


-V3 




(3,-2) 


3\/6 




(3,-3) 


V7 



Table Al. Coefficients for the quadratic "mixed" £ = 3 map 



Equation 


Term 


Coefficient 


Equation 


Term 


Coefficient 


4 


(13) 


-V42 


3 


(03) 


3V7 




(22) 






(12) 


%/T4 


2 


(02) 




1 


(01) 


%/l5 




(11) 


VTo 




(-12) 


-4V2 




(-13) 


-6\/3 




(-23) 


-Vso 





(00) 


3 










(-11) 


-1 










(-22) 


-7 










(-33) 


-3 









Table A2. Coefficients for the quadratic "mixed" £ = A map 
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Equation 


Term 


Coefficient 


Equation 


Term 


Coefficient 


4 


(04) 


14 


3 


(-14) 


7a/10 




(13) 


-7\/l0 




(03) 


-21 




(22) 


3\/70/2 




(12) 


V70 


2 


(-24) 


3V70 


1 


(-34) 


7\/T0 




(-13) 


-V70 




(-23) 


V70 




(02) 


— 11 




(-12) 


-6\/l0 




(11) 


3\/T0 




(01) 


9 





(00) 


9 










(-11) 


-9 










(-22) 


-11 










(-33) 


21 










(-44) 


14 









Table A3. Coefficients for tlie quadratic "pure" £ = 4 map 



Equation 


Term 


CoefRcient 


Equation 


Term 


Coefficient 


3 


(-333) 


-45 


2 


(-323) 


-45 




(-223) 


45 




(-213) 


10yi5 




(-113) 


-15 




(-222) 


20 




(-122) 


-5\/T5 




(-103) 


-15\/2 




(003) 







(-112) 


-35 




(012) 


15\/2 




(002) 


30 




(111) 






(Oil) 


^30 


1 


(-313) 


-15 





(-303) 







(-322) 


-5\/T5 




(-312) 


-15^2 




(-203) 


-15V2 




(-202) 


60 




(-212) 


35 




(-211) 


-5\/6 




(-1-13) 


-6VT5 




(-101) 


-36 




(-102) 


2^30 




(000) 


18 




(-111) 


-41 










(001) 


18 









Table A4. Coefficients for tlie cubic "pure" £ = 3 map C^^) 
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Equation 


Term 


Coefficient 


Equation 


Term 


Coefficient 


4 


(-444) 


30/7 


3 


(-434) 


30/7 




(-334) 


-30/7 




(-324) 


-5\/7/2 




(-224) 


-5/7 




(-333) 


5/56 




(-233) 


5V7/4 




(-214) 


15/4 




(-114) 


45/14 




(-223) 


135/56 




(-123) 


-15/4 




(-104) 


—v/TO/4 




(004) 


-13/7 




(-113) 


-65/56 




(013) 


\/T0/4 




(-122) 


-45V7/56 




(022) 


3^5/14 




(003) 


-17/28 




(112) 


-5\/7/14 




(012) 


23/4y^5/14 










(111) 


-15/28V7 


2 


(-424) 


-5/7 


1 


(-414) 


-45/14 




(-433) 


5\/7/4 




(-423) 


15/4 




(-314) 


-15/4 




(-304) 


-\/T0/4 




(-323) 


-135/56 




(-313) 


-65/56 




(-204) 


3^10/7 




(-322) 


-45\/7/56 




(-213) 


45\/7/28 




(-2-14) 


5V7/7 




(-222) 







(-203) 


23/4y^5/14 




(-1-14) 


-10\/7/28 




(-212) 


15/8 




(-103) 


-23/4^5/14 




(-1-13) 


-45\/7/28 




(-112) 


-15/8 




(-102) 


3\/T0/4 




(002) 


7/2 




(-111) 


-25/8 




(Oil) 


-3a/T0/8 




(001) 


1 





(-404) 


-26/7 








(-413) 


\/lO/4 








(-422) 


3a/70/14 








(-3-14) 


VTo/A 








(-303) 


17/14 








(-312) 


-23/4a/5/14 








(-2-24) 


3a/5/14 








(-2-13) 


-23/4^5/14 








(-202) 


7 








(-211) 


3\/T0/8 








(-1-12) 


3^10/8 








(-101) 


-2 








(000) 


1 







Table A5. Coefficients for tlie cubic "pure" i = 4 map 
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